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ABSTRACT 

Aims. UVES spectra of the very young ( ~ 10^ years) peculiar B-type star HR6000 were analyzed in the near-UV and visual spectral 
regions (3050-9460 A) with the aim to extend to other spectral ranges the study made previously in the UV using lUE spectra. 
Methods. Stellar parameters reff= 12850 K, log g=4. 10, and ^=Okm s"', as determined from H^, Hy, Balmer profiles and from the 
Fei, Fell ionization equilibrium, were used to compute an individual abundances ATLAS 12 model. We identified spectral peculiarities 
and obtained final stellar abundances by comparing observed and computed equivalent widths and line profiles. 
Results. The adopted model fails to reproduce the (b-y) and c color indices. The spectral analysis has revealed: the presence of 
emission lines for Mnn, Crn, and Pen; isotopic anomalies for Hg , Ca ; the presence of interstellar lines of Nai at AA 3302.3, 3302.9, 
5890, 5896 A, and of Ki at 7665, 7699 A; the presence of a huge quantity of unidentified lines, which we presume to be mostly 
due to Fen transitions owing to the large Fe overabundance amounting to [-1-0.7]. The main chemical peculiarities are an extreme 
overabundance of Xe, followed by those of Hg, P, Y, Mn, Fe, Be, and Ti. The most underabundant element is Si, followed by C, N, 
Al, S, Mg, V, Sr, Co, CI, Sc, and Ni. The silicon underabundance [-2.9] is the lowest value for Si ever observed in any HgMn star. 
The observed lines of Hei can not be reproduced by a single value of the He abundance, but they require values ranging from [-0.8] 
to [-1.6]. Furthermore, when the observed and computed wings of Hei lines are fitted, the observed line cores are much weaker than 
the computed ones. From the present analysis we infer the presence of vertical abundance stratification for He, Mn, and possibly also 
P. 
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1. Introduction 

HR6000 (HD 144667) is one of the most remarkable chemi- 
cally peculiar (CP) stars. It does not fit any of the CP sub- 
classes, but it seems to combine abundance anomalies from 
a variety of Bp sub-types. It forms with the star HR5999 
(HD 144668) the common proper motion visual binary system 
A199 or Dunlop 199 (Bessell & Eggen, 1972). The angular sep- 
aration between the HR 6000, which is the brighter more mas- 
sive component, and the secondary component, the well-known 
Herbig Ae star HR 5999, is about 45". 

The main interest to study the chemical composition of 
HR6000 comes from the generally poor understanding of the 
occurrence of abundance anomalies in such a young object with 
an estimated age of the order of 10^ years. In fact, the A 199 dou- 
ble system is located close to the center of the Lupus 3 molec- 
ular cloud which is populated by numerous TTauri stars. This 
has led to an assumption that this system has the same age as 
the cloud, which is estimated to be (9.1+2.1)xl0^ years (James 
et al., 2006). However, after the Hipparcos mission, the mem- 
bership of HR 5999 and HR 6000 to the cloud became rather 
questionable due to the uncertainties in the distance determina- 
tion for the Lupus cloud. The distance of HR 5999 and HR 6000 
measured by Hipparcos are 208+38 pc and 240+48 pc, respec- 
tively, while the Lupus cloud distance is 150+lOpc according 
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to Crawford (2000). Comeron et al. (2003) assigned a distance 
of about 200 pc to the Lupus cloud, but this determination was 
made by assuming a priori that A199 belongs to the cloud. 

A ROSAT survey of Herbig Ae/Be stars presented by 
Zinnecker & Preibisch (1994) has detected strong X-ray emis- 
sion coming from the direction of HR 6000. To explain the X-ray 
origin, a possible T Tauri companion for HR 6000 was postulated 
by van den Ancker et al. (1996) on the basis of an infrared ex- 
cess in the energy distribution relative to predictions made for 
an effective temperature of re(f= 14000 K. Siebenmorgen et al. 
(2000) fitted very short spectral regions of observed flux in the 
mid IR to a black body of 13000K and noticed that the observed 
spectrum was featureless. 

Neither spectrum nor radial velocity variability was found 
for HR6000 by Andersen & Jaschek (1984), who studied 
optical spectra taken at different epochs. However, van den 
Ancker et al. (1996) discovered long-term variations in the u, 
V, and b Stromgren magnitudes with approximate amplitude of 
0^03+07^01 in u, 0^02±0?'01 in v, and 0^01 in b. No variations 
in y were found. Kurtz & Marang (1995) discovered variations 
of 0.008 mag in V with a period near to 2 d and suggested that 
this period could be possibly caused by rotational variation in a 
spotted magnetic star This would imply a pole on orientation for 
HR6000 as they measured v sin /<5 km s ^ 

Catanzaro et al. (2004) were the first ones who provided 
abundances from the optical range. Previously, optical spec- 
tra were studied only by Andersen & Jaschek (1984) and 
Andersen et al. (1984) who identified the spectra from 3323 A 
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to 5317 A and discussed the abundances on the basis of the 
line intensities. It is remarkable that the stellar parameters 
derived by Catanzaro et al. (2004) from the Balmer profiles 
(reff=12950+50 K, log ^=4.05+0.01) are lOOOK lower than 
those previously adopted and deduced from the photometry 
by Castelli et al. (1985) (reif=14000K, log g=4.0), by Smith 
(1997) (reft = 13990K, log g=4.29) and by van den Ancker et al. 
(1996) (reff=14000K, log g=4.3). 

In this paper we examine the whole optical spectrum of 
HR 6000 from 3050 A to 9460 A with the aim to extend the anal- 
ysis performed on lUE spectra by Castelli et al. (1985) to the vis- 
ible region and to investigate the possible contamination of the 
stellar spectrum by a close T Tau companion. We re-examine the 
parameter determination, the line identification, and the abun- 
dances. We also searched for the presence of peculiarities like 
emission lines and isotopic anomalies which were recently dis- 
covered in HgMn stars. Owing to the availability of spectra taken 
at different epochs we also investigated possible spectral vari- 
abilities. 

The comparison of the observed and computed spectra as 
described in this paper is available at the web-address given in 
the footnot^E 



2. Observations and Radial Velocity 

Spectra of HR 6000 were recorded at ESO on 2005 September 
19 (JD 2453632.53881) and on 2006 March 25 (JD 
2453819.90226). In both epochs we used the UVES DICl and 
DIC2 standard settings to cover the spectral range from 3050 A 
to 9500 A. The slit width was set to 0'.'3 for the red arm, corre- 
sponding to a resolving power of /l/A/l ^ 1 10000. For the blue 
arm, we used a slit width of 0'.'4 to achieve a resolving power of 
~ 80000. The spectra were reduced by the UVES pipeline Data 
Reduction Software (version 2.5; Ballester et al. 2000) and us- 
ing standard IRAF routines. The signal-to-noise ratio (S/N) of 
the UVES spectra is very high, ranging from 150 in the near UV 
(3300A) to about 350 and 300 at 5000A and 7100 A, respec- 
tively. 

There are two gaps in the observed range at AA 4523^769 A 
and 7536-7660 A, which are caused by the physical gap between 
the two detector chips of the red CCD mosaic. In addition, the 
reduction in the spectral range 8076-8093 A is affected by the 
presence of a bad column on the MIT CCD. 

The heliocentric radial velocities derived from spectra ob- 
served on September 2005 and on March 2006 are 2.6 km sec ' 
and 1.4 km sec"', respectively. These values are consistent with 
the mean heliocentric radial velocity +2.6+1.8 kms"' of the 
Lupus cloud (James et al., 2006), so that they strenghten the as- 
sumption of the membership of HR6000 to the Lupus cloud. 
However, Catanzaro et al. (2004) measured 0.67+0.38 km s"', 
while Andersen & Jaschek (1984) determined -1.5±0.5kms"'. 
In their study, the values obtained from different spectra range 
from +0.6+0.8 km s"' to -3.3+0.8 km s"'. Such differences in 
measured radial velocities by different authors indicate that 
small radial velocity variations in HR 6000 are possibly real. 

The spectra were normalized to the continuum by using a 
modified version of the interactive code NORMA written by 
Bonifacio (1989). In the original procedure smoothed interpo- 
lated curves through selected continuum points are obtained by 
means of Hermite spline functions (subroutine INTEP by Hill 
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(1982)). We replaced the spline functions with a linear interpo- 
lation because in some other previous stellar analyses performed 
by one of us (EC.) it has given a continuum more close to that 
we would have drawn by hand. However, we checked that in 
the case of HR 6000 the continua from the two different inter- 
polations do not differ in a significant way. When needed, the 
continuum was adjusted over 6 A intervals with the help of the 
synthetic spectrum. The continuous level just longward of the 
Balmer and Paschen discontinuities is highly uncertain because 
it is difficult to know where to put it. Also the continuum drawn 
over the broad hydrogen lines is uncertain owing to the not fully 
corrected distorsions of the echelle spectra. 

3. Stellar parameters 

Table 1 summarizes the stellar parameters as determined from 
Stromgren, UBV, and near-infrared IJ photometry. The observed 
Stromgren and UBV indices were taken from the Hauck & 
MermilHod (1998) Catalogufl while (V-I,.)=-0.05 was taken 
from the Hipparcos Catalogue and (V-J)=-0.032 was obtained 
from the J magnitude given by The et al. (1996). The param- 
eter determination is the result of an interpolation of the ob- 
served dereddened indices in the grid of synthetic indices based 
on the NEWODF model atmospheres (Castelli & Kurucz, 2004) 
computed for metallicity [M/H]=0.0 and microturbulent veloc- 
ity ^=0.0 km s"'. The color grids are available at the web-address 
given in the footnotfl For the Stromgren and UBV photometry 
the adopted grids were uvbybeta and UBV, while for the (V-Ic) 
and (V-J) indices the adopted grid was UBVRIJHKL. 

Results from the Stromgren photometry are listed in the first 
panel of Table 1. For E(B-V)=0.06, the parameters from the 
Co and /3 indices are reff=l 3827+200 K, log ^=4.27+0.05. The 
reddening E(B-V)=0.06 was determined by Bessell & Eggen 
(1972) and later confirmed by The & Tjin A Djie (1978). An 
uncertainty in E(B-V) of 0.02 mag corresponds to a difference 
in Teff of only 15 K. The reddening E(B-V)=0.042 listed in the 
first row of Table 1 corresponds to E(b-y)=0.031 as yielded by 
the UBVYLIST code of Moon (1985) on the basis of standard 
reddening relations. The conversion E(b-y)=0.74E(B-V) from 
Crawford & Mendwewala (1976) was applied here. 

A microturbulent velocity of 2.0 km s"' lowers the tempera- 
ture by 100 K and does not affect the gravity, while an increase 
in metallicity to [M/H]=+0.5 decreases T^ff by about 50 K; a 
decrease in metallicity to [M/H]=-1.0 increases Teff by about 
+300 K. The corresponding differences in log g are of the order 
of 0.01 dex in all the cases. These values are only indicative of 
the metallicity effect in that the abundance pattern of HR 6000 is 
very different from any scaled solar one. 

The parameters from the UBV photometry and /3 index 
are listed in the second panel of Table 1. For E(B-V)=0.06, 
[M/H]=0.0, and ^=0.0 kms"' they are Teff^ 13973 +300 K, 
log ^=4.38+0.2 (Fig. 1). They are close to those obtained from 
the Stromgren photometry. This parameter determination, de- 
scribed in Castelli et al. (1985), is very sensitive to small red- 
dening changes. In fact, for E(B-V)=0.04 or E(B-V)=0.08 there 
is no intersection in the curves (U-B)o, (B-V)o and /3. For B- 
type stars, the (U-B), (B-V), and /3 indices do not depend on 
the metallicity in a significant way (Castelli, 1999; Castelli & 
Kurucz, 2006). 

The third panel of Table 1 shows that for a prefixed grav- 
ity log g=4.0, the observed (V-I)c and (V-J) indices are repro- 

" http://www.unige.ch/gcpd/gcpd.html 
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Fig. 1. Determination of the atmospheric parameters (Teffjog g) 
from (U-B)o, (B-V)o and fi. The circle represents the param- 
eters reff=13873K, log g=4.38 obtained from this method for 
[M/H]=O.Oand ^=O.Okms ^ The observed indices were dered- 
dened for E(B-V)=0.060 



duced by Teff equal to 10889 K and 1 15 17 K, respectively, when 
[M/H]=0.0, ^=O.Okms-', and E(B-V)=0.06. We used the red- 
dening relations E(V-I,)=1.25E(B-V) from Dean et al. (1978) 
and E(V-J)=2.23E(B-V) from Mathis (1990). The temperature 
was derived by interpolating the observed dereddened (V-Ic)o 
and (V-J)o indices in the Bessell et al. (1998) color grid com- 
puted for [M/H]=0.0 and ^=0.0 km s"' . 

Catanzaro et al. (2004) were the first ones who derived pa- 
rameters for HR 6000 from high resolution spectra, in particu- 
lar from Hy and Balmer profiles observed on FEROS spec- 
tra. They fixed Teff^l 2950+50 K and log §=4.05+0.01 from 
ATLAS9 models computed for [M/H]=-0.5 and zero helium 
abundance. This particular grid of models yields the same pa- 
rameters from both Balmer profiles, but we note that the param- 
eters do not change if the grid computed with [M/H]=0.0 and 
solar He abundance is used, provided that the different parame- 
ters from Hy ([12750K,3.95] and ([13150K, 4.13]) are av- 
eraged. As discussed in Sect. 4, parameters very close to those 
of Catanzaro et al. (2004) were derived by us from Balmer pro- 
files and Fe i, Fe ii equivalent widths measured on UVES spectra 
and from model atmospheres computed for the individual abun- 
dances of HR 6000. 

The last panel of Table 1 shows that the cooler 
model, as adopted by Catanzaro et al. (2004) (model 
T12950g405k2M05nohe), does not reproduce the observed 
Stromgren indices. In fact, the difference between the observed 
and computed (b-y)o indices amounts to 0.020 mag, or even 
more, depending on the adopted reddening, while the difference 
between the Cq indices is 0.122 mag. These discrepancies do not 
change in a remarkable way for the ATLAS 12 model (model 
T12850g41k0atl2) adopted by us to predict the observed 
spectrum (see Sect. 4). On the other hand, the hotter model 
[13900,4.32] fails to reproduce all the Balmer profiles. For solar 
abundances it provides good agreement between observations 
and computations for Hq., but the agreement decreases from H^ 



to Ha in the sense that the computed wings become too broad 
as compared with the observed ones; for individual abundances 
the computed inner wings of all the profiles are narrower than 
the observed ones, in particular for Hq, whose profile seems also 
to be a little bit asymmetric with a more extended red wing. All 
the comparisons between the observed and computed Balmer 
profiles are available at the website of footnote 1. Fig. 2 shows 
the case of H/j. 

In conclusion, there is a discrepancy of about lOOOK be- 
tween parameters derived from UBV and Stromgren photome- 
try and parameters derived from the line spectrum. Furthermore, 
the difference in T^s from the near-IR photometry is so large that 
spectrophotometric observations are needed in order to confirm 
these results. Unfortunately, the only spectrophotometric data 
available for HR6000 are the lUE spectra. For E(B-V)=0.06, 
the best fit of the UV energy distribution from the lUE images 
SWP14849 and LWR11431 to the grid of ATLAS9 energy dis- 
tributions computed for [M/H]=0.0 and ^=0.0 km s ' is achieved 
with Teff = 1 3300 K, log §=4. 1 . More details about these lUE im- 
ages can be found in The et al. (1996). 

4. The ATLAS12 model and abundance 
determination 

The opacity distribution function ATLAS9 model with pa- 
rameters derived from the UBV and Stromgren photometry 
(reff = 1 3900 K, log §=4.32, [M/H] =0.0, ^=0.0 km s" ' ) was used 
to derive preliminary abundances. Abundances for iron and 
phosphorus were obtained from equivalent widths, while those 
for all the other elements were obtained with the synthetic 
spectrum method. Mostly the lines listed in Castelli & Hubrig 
(2004b) (electronic Appendix) were examined. Synthetic spec- 
tra were computed with the SYNTHE code (Kurucz, 1993). The 
adopted line lists are available on the web-address given in the 
footnot^. The synthetic spectra were broadened for a rotational 
velocity vsin/=1.0kms"' and an instrumental resolving power 
of 80000 for /}<4520 A and of 1 10000 for /}>4765 A (Sect. 2). 
The adopted v sin / was that reproducing the observed profiles at 
best, in particular Mg ii 448 1 A. It is in close agreement with the 
value of 0.0 km s"' derived by Catanzaro et al. (2004). 

With this model, the average abundance derived from the 
equivalent widths of 11 Fei selected lines was equal, within 
the errors limits, to that obtained from the equivalent widths 
of 34 selected Fen lines, i.e. log(N(Fei)/Nto,)=-3.65+0.06 
and log(N(Fen)/Nfo,)=-3.69+0.15. At the same time, for 
^=0.0 km s"' there was no trend of the abundances versus the 
equivalent widths of the Feii lines. The table with the measured 
equivalent widths and corresponding abundances is available at 
the website of footnote 1 . 

The abundances based on the ATLAS9 model were used as 
input data for an opacity sampling ATLAS 12 model (Kurucz 
1997; Castelli & Kurucz, 1994) computed for the same pa- 
rameters of the ATLAS9 model. Due to the extremely peculiar 
abundances of HR6000, the 'Y-tjioss structure of the ATLAS 12 
model is rather different from that of the ATLAS9 model, in par- 
ticular in layers upper than log(TRosj)< -0.5. The comparison 
of the T-tross relations from ATLAS9 and ATLAS 12 is simi- 
lar to that shown in Castelli & Hubrig (2004b). As a conse- 
quence, in spite of the achieved better agreement between ob- 
served and computed Balmer profiles, the Fe i-Fe ii ionization 
equilibrium worsens. The average abundances for Fei and Feii 
became -3.42+0.06 dex and -3.61+0. 15 dex, respectively. 

http://wwwuser.oat.ts.astro.it/castelli/linelists.htmI 
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Table 1. Model parameters for HR 6000 from photometry by assuming [M/H]=0.0 and ^=0.0 km s"^ . 



Observed indices Dereddened indices Parameters 


(b-yj nil C| 13 


l-(b-y) (b-y)o ni,, c,, 


TeiiCK) log g l^(B-V) 


-0.030 0.116 0.512 2.750 
±0.003 ±0.003 ±0.003 ±0.013 


0.031 -0.061 0.126 0.506 
0.044 -0.074 0.130 0.503 


13799 4.27 0.042 
± 150 ±0.05 
13827 4.27 0.060 
±200 ± 0.05 


V U-B B-V 


E(B-V) Vo (U-B)o (B-V)o 


reff(K) log g 


6.647 -0.434 -0.073 
±0.013 ± 0.004 ± 0.004 


0.06 6.461 -0.477 -0.133 


13973 4.38 
±300 ±0.2 


(V-I), (V-J) 1 E(V-1,) E(V-J) (V-l,)o (V-J)o 1 r,ff(V-J) reff(V-I,) log g 


-0.05 -0.032 1 0.075 0.134 -0.125 -0.166 | 11517 10889 4.0 (prefixed) 
±0.00 


model 


O'~y)co/c ^calc 




T12950G405k2M05nohe 
T12850g41k0atl2 


-0.041 0.629 
-0.049 0.596 






i?oo HOO 160* laflo ?o83 nm !Mt) jmb 3000 

Fig. 3. Comparison of lUE data (dotted line) with the flux 
computed from the ATLAS 12 model with individual stellar 
abundances (full line). Model parameters are refF= 12850 K, 
log ^=4.1, and ^=0.0kms The lUE flux was dereddened for 
E(B-V)=0.06. 



A search for an appropriate ATLAS 12 model suitable to re- 
produce both the Balmer profiles and the Fei-Fen ionization 
equilibrium has finally yielded, after several trials and errors, the 
parameters reff= 12850 K, log g-A.\, in good agreement with 
the values adopted by Catanzaro et al. (2004). The Fei and Fen 
abundances were -3.89±0.06dex and -3.88±0.16dex, respec- 
tively. Fig. 3 shows that this model also reproduces, for E(B- 
V)=0.06, the slope of the UV energy distribution given by the 
lUE images SWP14849 and LWR11431. However, a close in- 
spection of Fig. 3 reveals the presence of some strong absorp- 
tions in the short-wavelength region which are neither predicted 
by the models nor detected in the examined lUE high-resolution 
spectra (Sect. 5.1). As already discussed in the previous sec- 
tion, the final ATLAS 12 model fails in reproducing the observed 
Stromgren indices. The conclusion is that we were not able to 
fix a model well suited to reproduce the atmosphere of HR 6000 
as deduced from the available observations. Table 2 summarizes 
the parameter values from the different determinations based on 
flie ATLAS9 (AT9) and ATLAS12 (AT12) models. 



Table 2. Sunnmary of the stellar parameters for HR 6000 from 
different determinations. 



Method 




log g 


Models 


UBVjS 


13973±300 


4.4±0.2 


AT9 




13813±200 


4.3±0.05 


AT9 


lUE flux 


13300±100 


4.1±0.1 


AT9 


Balmer profiles 


12850±50 


4.1±0.01 


AT12 


Fe i/Fe 11 


12850 


4.1 


AT12 


VIJ 


11203±300 


4.0 (prefixed) 


AT9 



5. The resulting abundances 

The abundances used for computing the whole synthetic spec- 
trum of HR6000 from 3050 A to 9460 A are listed in col- 
umn 2 of Table 3. The synthetic spectrum overimposed on the 
observed spectrum is available at the website of footnote 1. The 
abundances are based on the ATLAS 12 model with parameters 
refF=12850K, log ^=4.1, ^=O.Okms-i, and vsinj=1.0kms-i. 
They are shortly discussed below. 

The main results are the large underabundance of Si ([-2.9]) 
followed by that of C ([-1.8]). Other underabundant elements 
are He (<[-0.85]), N < ([-1.7]), O ([-0.5]), Ne([-0.9], Mg 
([-1.2]), Al ([-1.7]), S ([-1.5]), CI (<[-0.7]), Ca ([-0.2]), Sc 
([-0.6]), V ([-1.1], Co < ([-0.8]), Ni ([-0.4]), and Sr ([-1.0]). 

On the other hand, the most overabundant element is Xe 
([-1-4.5]), foUowed by Hg ([-1-2.7]), and P ([+2.3]:). Other over- 
abundant elements are Be ([-1-0.7]), Ti ([+0.25]), Mn ([+0.95]), 
Fe ([+0.7]), Y ([+1.0]), and possibly Ba([+1.0]). 

Below we discuss the individual abundances: 

Helium: From the analysis of the Hei lines 4713 A and 
5875 A, Catanzaro et al. (2004) derive an underabundance of 
[-1.45], i.e. N(He)/N„„=:0.0032. We extended their analysis 
to all the He I lines previously analyzed in the HgMn star 
HD 175640 by CastelH & Hubrig (2004b), except for Hei 
4713 A which is located in a gap of the observed spectral 
regions. While all the Hei lines were fairly well predicted 
in HD 175640 by the same abundance, this is not a case in 
HR6000. In fact, the abundance N(He)/Ntoj from the different 
lines ranges from 0.0125 to 0.002, i.e from [-0.80] to [-1.60]. 
The higher value reproduces the whole profile of Hei 3867.4 A 
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Fig, 2. The observed (thick hne) is compared with profiles computed with different models (thin lines). In the upper and mid- 
dle panels the ATLAS9 model and the ATLAS12 model with parameters given by the photometry (reff=13900K, log g=4.32, 
^=0.0 km s ') were used. The lower panel shows the profile computed with the final adopted ATLAS 12 model with parameters 
Teff = 1 2850 K, log g=4. 1 , and ^=0.0 km s" ' . 



and only the wings of He i lines at 4026. 1 A and 4387.9 A, whose 
computed cores appear stronger than the observed ones (Fig. 4). 
The same discrepancy occurs also for the lines at AA 4471.4, 
4921.9, 5015.7 A, whose wings are fitted by abundances equal 
to 0.005 and 0.006, while the cores would require a lower value. 
The two weak lines at AA 5875.6 and 6678 A are rather well 
predicted for N(He )/N,or equal to 0.002 and 0.003, respectively. 
Finally, the line at 7065.2 A is partially blended both with tel- 
luric lines and with Fen 7065.163 A. It can not be reproduced 
whichever is the abundance. The abundance from the remaining 
He I lines, which are very weak, depends on the assumed posi- 
tion for the continuum level. The comparison between observed 
and computed He i profiles is available at the website of footnote 
1. 

Profiles of He I like those observed in HR6000 were ob- 
served also in other HgMn stars as for instance HR7361, 
HR7664 (Dworetsky, 2004) or Feige 86 (Bonifacio et al., 1995). 
This kind of anomalous profiles was interpreted as evidence of 
helium stratification (Dworetsky, 2004; Bohlender, 2005). 

No contribution of ■'He to A 6678 A was observed. 

Beryllium: The abundance was derived from the lines at 
3130.42 and 3 13 1.06 A, although the second one is heavy 
blended. 

CNO: While Catanzaro et al. (2004) did not observe Cii 
lines, we possibly identified very weak C ii lines at AA 3918.968, 
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Fig. 4. Comparison of observed (thick line) and com- 
puted (thin line) Hei 4026 A profiles. The abundance 
log(N(He)/N,o,)=-1.90 which fits the wings yields a too strong 
computed core. 
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Table 3. Final abundances log(Neiem/Ntot) from the ATLAS 12 
model with parameters reff=12850K, log g=4.l, ^=O.Okms"'. 
Solar abundances are from Grevesse & Sauval (1998). 



Elem 


UVES 


lUE(old) 


lUE(new) 


Sun 


He 


-1.90 






-1.05 


Be 


-9.90 


<-9.5 


-9.90 


-10.64 


B 




-10.1 


-10.29 


-9.49 


C 


-5.32: 


>-5.4 


-5.82 


-3.52 


N 


< -5.82 


>-6.5 


-6.50 


-4.12 





-3.71 


-3.3 


-3.70 


-3.21 


Ne 


-4.86 






-3.96 


Na 


-5.71 






-5.71 


Mg 


-5.66 


-5.6 


-5.66 


-4.46 


Al 


-7.30 


-7.7 


-7.30 


-5.57 


Si 


-7.40 


<-5.7 


-7.40 


-4.49 


P 


-4.30 


-4.4—5.0 


-5.00 


-6.59 


S 


-6.20 


-6.0—5.5 


-6.20 


-4.71 


CI 


<-7.24 


-4.3 


-7.24 


-6.54 


Ca 


-5.88 


-5.2 


-5.88 


-5.68 


Sc 


-9.50 


-11.0 


-9.50 


-8.87 


Ti 


-6.77 


>-6.5 


-6.77 


-7.02 


V 


-9.14 


-6.7—8.0 


-9.14 


-8.04 


Cr 


-6.37 


-6.2 


-6.37 


-6.37 


Mn 


-5.60 


-5.3 


-5.60 


-6.65 


Fe 


-3.85 


-3.8 


-3.85 


-4.54 


Co 


<-7.92 


>-9.3 


-7.92 


-7.12 


Ni 


-6.19 


-7.0 


-6.19 


-5.79 


Cu 


< -7.83 


>-9.3 


-7.83 


-7.83 


Zn 




-8.3 


-7.44 


-7.44 


Ga 


<-9.16 


-8.4 


-9.16 


-9.16 


Sr 


-10.07 






-9.07 


Y 


-8.80 






-9.80 


Xe 


-5.40 






-9.87 


Ba 


-8.91? 






-9.91 


Hg 


-8.20 




-8.20 


-10.91 



4267.01, and 4267.26 A. No nitrogen hues were observed. We 
deduced a carbon underabundance of [-1.80] and an upper limit 
for nitrogen of [-3.40]. Oxygen is underabundant by [-0.50]. We 
examined all the Oi lines listed in Castelli & Hubrig (2004b). 
The Oi infrared lines at AA 111 A, 8446, and 9260 A are inade- 
quately computed under the LTE assumption. 

Neon: The abundance was fixed from Ne i 7032.413 A. 

Sodium: There is a remarkable inconsistency between the 
abundance of -3.31 dex from the two lines Nai at 3302.368 A 
and 3302.978 A and the almost solar abundance of -5.57 dex of 
the Na I lines at AA 5688.2, 8183.2 and 8194.8 A. As already de- 
scribed by Catanzaro et al. (2004) the two stellar Nai D lines 
at 5889.951 A and 5895.924 A are completely hidden by the in- 
terstellar lines. The most likely explanation is that the lines at 
3302 A are also of interstellar origin. They have also been ob- 
served in the spectrum of HD 175640 as weak red-shifted ab- 
sorptions with the same wavelength displacement of -1-0.08 A as 
the interstellar D-lines. 

Magnesium: The estimated underabundance is [-1.20]. 
Magnesium abundance was mostly derived from Mg n 448 1 A. 
The Mg II doublet at 6545.942, 6545.994 A is much broader than 
predicted and is also red-shifted. The most plausible explanation 
is the lack of Stark broadening data in the computations. 

Aluminium: The estimated underabundance is [-1.73] from 
the Ahi fines at AA 7042.06 and 7056.60 A. No otfier fines, as 
Al n at AA 3944.009 and 3961.523 A, were observed. 



Silicon: The estimated underabundance is [-2.90]. Only 
weak Sin fines at AA 3853.665, 3856.018, and 3862.595 A were 
observed. 

Phosphorus: The overabundance of phosphorus can be in- 
ferred from most of the Pn and Pm lines, but the real value 
is very difficult to be deduced because the different lines give 
rather different abundances. We measured the equivalent widths 
for a sample of 22 Pn and 4 Pin lines. The average abundance 
from P n is log(N(P)/Nto,)=-4.36+0.28 for ^=0.0 km s"' , but the 
trend of the abundances with the equivalent widths is very large. 
It disappears for ^=6.0 km s"^ The corresponding abundance 
of -4.70+0. 12 produces, for f =6.0km s"', extremely broad and 
shallow profiles completely different from the observe ones. 

The abundance for ^=0.0 km s ' inferred from the Pin lines 
is log(N(P)/Nto,)=-4.57±0.27. Tfierefore, tfie Pn abundance for 
^=0.0 km s"^ is 0.2 dex larger than that obtained from Pm, but 
lies withing the large mean square error of the averages. The 
large uncertainty in the phosphorous abundance could be related 
to NLTE effects which have been neglected here. Another pos- 
sibifity could be tfie presence of pfiospfiorus vertical abundance 
stratification. 

Tfie adopted P n and P in lines, togetfier witfi tfie atomic data 
and tfie measured equivalent widtfis, are fisted on tfie website of 
footnote 1. 

Sulfur. The sulfur underabundance of [-1.5] was derived 
from the same lines used by Castelli & Hubrig (2004b). 

Chlorine: Tfie absence of observed cfilorine lines, like Cln 
4794.54 A, impfies an abundance less tfian -7.24 dex. 

Potassium: Potassium can be predicted by adopting solar 
abundances. However, tfie sligfit blue-sfiift of -0.05 A of tfie 
fines Ki AA 7664.911, 7698.974 A, tfie only ones observed in 
tfie spectrum, may indicate an interstellar origin for tfiem. 

Calcium: All the Can lines listed in Castelli & Hubrig 
(2004b) can be reproduced by an underabundance of [-0.2], 
except for the fines of mult. 13 (8201.720A and 8248.796A) 
wfiicfi are predicted too strong. The blue wing of the K line 
at 3933. 664A presents an additional broad weak absorption. It 
is further discussed in Sect. 6.3. The lines of tfie Can infrared 
triplet are red-sfiifted by 0. 14 A. Tfiey are discussed in more de- 
tail in Sect. 6.1. 

Tfie Ca I line at 4226.728 A is very weak. It is reproduced by 
solar abundance. A possible interstellar contribution can not be 
excluded. 

Scandium: Tfie underabundance of [-0.6] was derived from 
tfie Sc n lines listed in Castelfi & Hubrig (2004b). 

Titanium: The overabundance of [+0.25] was derived from 
the Tin lines fisted in Castelfi & Hubrig (2004b). 

Vanadium: Tfie underabundance of [-1.1] was deduced from 
tfie lines of Vn mult. 1 at 3100 A. 

Chromium: Most of tfie observed Crn lines can be predicted 
by solar abundance. 

Manganese: Tfie abundance of -5.0dex from tfie lines of 
Mnn mult.3 at 3440-3498 A is larger by 0.6 dex than that ob- 
tained from most of the Mn n lines lying longward of the Balmer 
discontinuity. This result indicates the possible occurrence of 
vertical stratification for tfie abundance of tfiis ion. Tfie abun- 
dance of -5.6 dex listed in Table 3 is the one obtained from the 
region longward of the Balmer discontinuity. Its value is the 
same as that presented by Catanzaro et al. (2004). 

Iron: The abundance was derived from the equivalent widths 
of 1 1 Fe I and 34 Fe ii lines. The lines togetfier witfi tfie atomic 
data and the measured equivalent widths are available on the 
website of footnote 1 . The overabundance of [+0.7] gives rise 
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in the spectrum to a very large number of Fe i, Fe iii, and espe- 
cially Fe II lines (see sub-section 6.4). 

Cobalt: No Co ii lines were observed in the spectrum, in par- 
ticular the line at 3501.717 A is absent. Therefore only an upper 
limit of [-0.8] was fixed. 

Nickel: Most Niii lines indicate an underabundance of 
[-0.4]. 

Copper: No lines of Cu ii were observed, in particular those 
at AA 4909.734, 4917.892, and 4931.698 A. Because for solar 
abundance there are no predicted lines we assumed solar abun- 
dances as upper limit for this element. 

Gallium: No gallium lines were observed. If a feature at 
A 5416. 32 A is identified as Gaii it would be reproduced by 
log(N(Ga )/N,o,)=-6.0. This abundance gives rise to other strong 
computed Ga ii lines with no observed counterparts. 

Strontium: No Srii lines were observed. Only a very weak 
absorption is detectable at the position of Srii 4077.71 A. It is 
reproduced by assuming an underabundance of [-1.0]. 

Yttrium: The overabundance of [-1-1.0] was derived mostly 
from the Y ii line at 4900. 12 A. 

Zirconium: No Zr lines were observed in the spectrum. They 
are not predicted by assuming solar Zr abundance. 

Xenon: Xenon is the most overabundant element ([-1-4.5]). 
The abundance was derived from the Xeii lines listed in 
Castelli & Hubrig (2004b) having log gf available from the 
NIST database. They were taken from Wiese & Martin (1980) 
(WM80). Furthermore, we added in the line list some lines with 
log g/'s from Ryabchikova & Smirnov (1989). In addition, only 
for identification purposes, all the lines with no log gf in the 
NIST database but with intensity > 100 were included in the 
line lists. A guessed log g f was assigned to them. 

We noticed that the position of several Xe ii lines is close to 
that of unidentified stellar lines with the wavelength blue-shifted 
by 0.05 to 0. 1 A from the Xe ii laboratory wavelength. Examples 
are the hnes at 4180. 10, 4208.48, 4209.47, 4238.25, 4245.38, 
4330.52, 4393.20, 4448.13, 4462.19, 5667.56A. We have ob- 
served the same occurrence also in the two Xe rich HgMn 
stars Feige 86 and 46 Aql (Castelli & Hubrig, in preparation). 
If the Xe II identification is correct, the wavelength shift could 
be due to some isotopic anomaly. In fact, no isotopic compo- 
sition was considered in our computations owing to the lack 
of isotopic wavelengths for Xeii. There are data only for A 
6051. 150A (Alvarez et al., 1979). For this line the isotopic 
wavelengths range from 605 1 . 148 A to 605 1 . 152 A. Because the 
stellar line is blue-shifted by 0.03 A from the laboratory wave- 
length 605 1.1 50 A (Hansen & Persson, 1987) some isotopic 
anomaly can not be invoked to explain the observed shift. 

Barium: We were not able to draw any conclusive results for 
barium. There is a weak absorption at 4933.93 A, but if we iden- 
tify it as Ba ii 4934.076 A, this barium line would be blue-shifted 
by 0.1 5 A. The wavelength in the spectrum would correspond 
to the Ritz wavelength 4933. 970A rather than to the observed 
wavelength 4934.077 A, as they are listed in the NIST database. 
Any isotopic anomaly has to be excluded because the shortest 
isotopic wavelength is that of the hyperfine component of '^^Ba 
at 4934.054 A (Mc William, 1998). 

Mercury: Only Hgii at 3984 A is observable. It can be pre- 
dicted by assuming an overabundance of [-1-2.7] and a non- 
terrestrial isotopic composition. The Hg isotopic anomaly is dis- 
cussed in Sect. 6.1. 



5. 1 . Comparison with the abundances from the lUE spectra 

Column 3 of Table 3 lists the abundances that were derived from 
lUE spectra by Castelli et al. (1985). The analyzed images were 
SWP04588 and LWR03979. We revised these abundances by us- 
ing the same ATLAS 12 model adopted for the analysis of the 
UVES spectra and the same images that we downloaded from 
the MAST archiv^ Also updated line lists, as compared with 
those used by Castelli et al. (1985), were used. Column4 of 
Table 3 lists the revised abundances from the ultraviolet. 

The lUE spectra are very noisy and the resolving power is 
too low (R=~ 12000- 13000) to allow one to perform an accu- 
rate analysis. Almost all very numerous lines are blended either 
with known or with unknown components. In addition, the never 
adequately corrected distortions of the echelle orders make the 
drawing of the continuum a very problematic task. The large un- 
certainty in the continuum position affects both the shape and 
the strength of the spectral lines. For all these reasons the abun- 
dances that we derived from the lUE spectra have always been 
considered as estimates rather than definitive values. 

Table 3 shows that the difference IUE(new)-IUE(old) be- 
tween revised and Castelli et al. (1985) abundances is larger 
than 0.4dex for C, O, CI, Ca, Sc V, Ni, Zn, and Ga. The 
abundances from UVES spectra reproduce the lUE spectrum 
(within all the above outlined limits) for all the elements except 
for carbon and phosphorous, for which the differences UVES- 
lUE(new) are 4-0.5 dex and +0.7 dex respectively. While the 
analysis of the carbon abundance from the visible spectral region 
is based only on very weak blended lines and therefore is rather 
uncertain, the phosphorous discrepancy confirms the difficulty 
in fixing the phosphorus abundance from the different lines. 

The synthetic spectrum overimposed on the lUE spectra is 
available at the website of footnote 1 . It was computed with the 
abundances Usted in column 4 of Table 3. 



6. Line spectrum peculiarities 

In the following sub-sections we discuss isotopic anomalies ob- 
served for Hg and Ca and the emissions observed for a few lines 
of Mn II, Cr ii, and Fe ii. We also discuss the huge quantity of 
lines which could not be identified. 

6. 1 . Isotopic anomalies 

Isotopic anomalies in HR 6000 were observed for the Ca ii in- 
frared triplet and for Hgn at 3984 A. No ^He/"^He isotopic 
anomaly was detected in the UVES spectra. 

Can infrared triplet: The lines of the infrared triplet at 
^^8498.023, 8542.091, and 8662. 141 A are red-shifted by 
0. 14 A and appear broader and stronger than the predicted ones. 
The wavelength red-shift of the infrared triplet observed in sev- 
eral CP stars (Cowley et al., 2007) was interpreted by Castelli 
& Hubrig (2004a) as due to a Ca isotopic composition different 
from the terrestrial one. 

Hgii 3984: Hgii at 3984 A, although weak, was detected in 
the UVES spectra. The isotopic composition that best repro- 
duces the observed line for log(N(Hg )/N,o,)=-8.20 is that listed 
in column 3 of Table 4. For comparison, the terrestrial isotopic 
composition (Anders & Grevesse 1989) is given in column 5 of 
Table 4. In HR6000 the isotope 204 is the most abundant one. 
The profiles computed with the stellar and terrestrial isotopic 
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Table 4. Mercury isotopic (iso) and hyperfine (hfs) composition 
in % for Hg II at 3984 A. 



Isotope 


A 


stellar 


log gf 


terr. 


log sf 


hfs. 






% 


iso 


% 


iso 


% 


196 


3983.769 


0.15 


-2.81 


0.15 


-2.81 




198 


3983.838 


0.77 


-2.11 


9.97 


-1.00 




199a 


3983.838 


20.00 


-0.70 


16.87 


-0.77 


7.14 


199b 


3983.849 


20.00 


-0.70 


16.87 


-0.77 


9.71 


200 


3983.909 


6.87 


-1.16 


23.10 


-0.64 




201a 


3983.930 


5.00 


-1.30 


13.18 


-1.32 


4.80 


201b 


3983.941 


5.00 


-1.30 


13.18 


-1.08 


8.30 


202 


3983.990 


10.00 


-1.00 


29.86 


-0.52 




204 


3984.071 


57.20 


-0.24 


6.87 


-1.16 
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Fig. 5. The observed profile (thick line) of Hg II 3984 A is com- 
pared in the lower panel with the profile computed with a terres- 
trial isotopic composition (thin line) and in the upper panel with 
the profile computed with the isotopic mixture listed in column 3 
of table 4 (thin line). 



compositions are compared with the observed profile in Fig. 5. 
We note that the Hg ii line is blended with Fe i 3983.956 A. 

Hg II at 1942.3 A can be observed in the lUE spectra, but the 
isotopic shifts are too small to enable us to confirm the isotopic 
composition deduced from the line at 3984 A. 

6.2. Emission iines 

Emission lines were observed for Mnii, Cm, and Feii. We as- 
sumed that the studied lines are true emissions only if they were 



observed in both UVES spectra taken at two different epochs. 
Emission lines are listed in Table 5. 

Emission lines of Mn II: Some lines of Mnii multi- 
plet 13 have been observed either in emission {AA 6125.863, 
6126.218, 6126.510, 6131.923 A) or much weaker than pre- 
dicted {AA 6122.434, 6122.810, 6128.734 A). Some others lines 
{AA 6129.033, 6130.796, 6131.016A) are not observed in the 
spectrum taken on September 2005, while they appear as weak 
emissions in the spectrum of March 2006. This fact may sug- 
gests the presence of some variability in the emissions. Finally, 
the two very weak Mnii lines of mult. 13 at 6123. 160 A and 
6129.254 A are in absorption and agree rather well with the pre- 
dicted lines. 

Other weak Mnii emissions occur at 6446.337 A (Mnn, 
mult. 19) which is blended with Feii 6446.410 A, and at 
7219.968, 8495.229, 8565.819, 8769.175 A. An emission could 
also affect the line at 8695.208 A, which is observed much 
weaker than predicted, and the line at 8819. 584A which, al- 
though predicted, is not observed. We note that there is a rather 
good agreement for this last line between the observed and com- 
puted profiles in the HgMn star HD 175640 (Castelli & Hubrig, 
2004b), even considering that no hyperfine structure is available 
for computations. 

Emission lines of Cr ii: A few Cr ii lines arising from high ex- 
citation levels and having large log gfs appear clearly in emis- 
sion. These lines were also found in emission in HD 175640 
(Castelli & Hubrig, 2004b). 

Emission lines of Fe ii: There are no clear emissions for Fe n. 
The most probable ones are listed in Table 5. However, there is 
a set of observed lines much weaker than predicted. This dis- 
agreement could be explained either due to incorrect log g/ s' or 
to emissions filling the absorptions. They occur at AA 6145.61 1, 
6214.948, 6328.498, 6340.841, 6349.601, 6419.626, 6436.598, 
6451.094, 6479.607, 6480.807, 6518.774, 6524.726, 6576.153, 
6576.170A. 

The phenomenon of weak emission lines at optical wave- 
lengths for main-sequence B-type stars has been noted from 
observations of relatively few stars. Presently, explanations of 
this phenomenon have been put forward in the context of the 
non-LTE line formation (Sigut, 2001) and possible fluorescence 
mechanisms (Wahlgren & Hubrig, 2000). 

6.3. Unidentified lines 

There is a huge number of unidentified lines in the spectrum 
of HR 6000. Their list is available at the web-address given in 
footnote 1 . For instance, most impressive regions are AA 4404- 
4411 A and 5100-5300 A, which are overcrowded with uniden- 
tified absorptions. We estimate that most of them are due to 
Fell. In fact, a comparison with the spectrum of HD 175640 
(Castelli & Hubrig, 2004b) has shown unidentified lines at the 
same wavelengths also in this star, although less intense than 
in HR6000, in accordance with the different iron abundances 
(Fig. 6). Therefore the unidentified lines are not a special pecu- 
liar characteristic of HR 6000. 

In the Kurucz line list for Fe ii (Kurucz, 2005 there are 
numerous transitions between observed levels with wavelengths 
close to those of the unidentified lines, but with so low log gf 's 
that none corresponding absorption is predicted. The compari- 
son of log gf 's from Kurucz with the log gf 's from Raassen 
& Uylings (1998) has shown agreement in the two determina- 
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Mult. 


1 f 
log Sj 




Notes 
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6112.906 


Fell? 


— 


-3.753 


90300.625 


not predicted. 


6122.434 


Mnii 


13 


0.950 


82136.400 


absorption filled by emission? 


6122.810 


Mnii 


13 


0.084 


82136.400 


absorption filled by emission? 


6125.863 


Mnii 


13 


0.783 


82144.480 




6126.218 


Mnii 


13 


0.230 


82144.480 




6126.510 


Mnii 


13 


-0.791 


82144.480 




6128.734 


Mnii 


13 


0.588 


o 1 c 1 1 /' r\ 

82151.160 


absorption filled by emission? 


6129.033 


Mnii 


13 


0.208 


82151.160 


variable? 


6130.796 


Mnii 


13 


0.354 


82155.840 


variable? 


6131.016 


Mnii 


13 


0.053 


82155.840 


blend with Fe n in absorption 


6131.923 


Mnii 


13 


0.053 


82158.170 




6158.621 


Cm 




0.718 


OrVITyl C\C}C\ 

89174.080 




6182.340 


Crn 




0.452 


or\ o f\f\ 

89336.890 




6446.337 


Mnii 


19 


-0.663 


98423.300 


on the blue wing of Fe n 6446.410 A 


6526.302 


Crn 


— 


0.173 


89885.080 




6585.241 


Crn 




0.829 


90850.960 


red-shifted 


7174.368 


Mnn 


— 


0.397 


85960.460 




7175.815 


Mnn 


— 


-0.023 


85960.460 




7219.968 


Mnn 




0.768 


86057.440 




7222.639 


Mnn 




0.252 


86057.440 




8495.229 


Mnn 


— 


0.198 


85960.460 




8565.819 


Mnn 


— 


0.278 


86057.440 




8695.208 


Mnn 




0.577 


74560.010 


absorption filled by emission? 


8740.349 


Fe n 




0.002 


99068.450 




8758.410 
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8769.175 


Mnn 




0.491 


74560.010 




8819.584 


Mnn 


13 


0.347 


74560.010 


absorption filled by emission? 


8839.005 


Fen 








J78, unclassified" 


9020.127 


Fen 








J78, unclassified" 


9085.11 


? 










9208.4 


9 











"J78=Johansson(1978) 



tions. Therefore, either both determinations are off by several 
orders of magnitudes or other transitions at the same wavelength 
occur. This possibility seems partly justified when the Kurucz 
line list including lines from both observed and predicted lev- 
els is usecfl. In fact, numerous Fe ii lines due to transitions be- 
tween predicted levels have log g/'s large enough to produce 
observable absorptions. Unfortunately, the uncertainty in their 
wavelength may be as large as 1 A and also the coiTespond- 
ing log ^/'s may be affected by rather large errors. Table 6 is 
an example of tentative line identification for the unidentified 
lines in the interval 5130-5136 A. As it is shown in Fig. 6 and in 
the lower panel of Fig. 7, out of displayed 17 spectral lines only 
three lines can be identified as arising from observed levels. They 
are Fe n 5 1 3 1 .2 1 A, Fe n 5 1 32.669 A, and Fe i 5 1 33 .688 A. The 
middle panel of Fig. 7 shows that more lines would be identified 
if several log 's listed in col. 5 of Table 6 are replaced by those 
given in col. 7. They were derived by fitting the computed lines 
to the observed lines for the stellar iron abundance of -3.85 dex. 
Upper panel of Fig. 7 shows that the number of identified lines 
increases if the Fe ii lines from predicted levels listed in cols. 8 
and 9 of Table 6 are added. The intensity of the lines from pre- 
dicted levels is usually lower than that of the observed lines. 

There is a remarkable coincidence between several Fell 
lines measured by Johansson (1978) (Table 6, col. 3) and the 
unidentified lines in our spectra. However, they are not predicted 



^ http://kurucz.harvard.edu/atoms/gf2601/gf2601.lines0600 



by the synthetic spectrum either because they are missing in the 
line lists or owing to their too low assigned log gf. 



7. The TTauri companion 

Van den Ancker et al. (1996) suggested the presence of a TTauri 
companion for HR6000 in order to explain the strong X-ray 
emission observed by Zinnecker & Preibisch (1994) and the in- 
frared excess that they inferred from the comparison of the ob- 
served energy distribution with that computed for the best fit- 
ting model having parameters retf=14000K, log §=4.3. They 
assumed a temperature of 3500 K for the T Tauri central star, an 
accretion rate of 2x10"^ Moyr"', and a V magnitude of about 
13.5, so that the companion would be about 6 mag fainter than 
HR 6000 and can not contribute significantly to the observed vi- 
sual spectrum. 

There are no clear spectroscopic traces of a second star in 
the spectrum of HR 6000, except for a broad weak absorption at 
XX (flQl-(flQ% A which, although of different intensity in the two 
spectra observed on September 2005 and March 2006, exhibits 
a very similar shape. The synthetic spectrum of HR 6000 con- 
sidered as a single star is flat in this region except for two weak 
Fe II lines at 6707.232 A and 6707.539 A which do not cover the 
whole observed absorption (Fig. 8, lower panel). We could argue 
that the feature is Li 1 6707.7 which is a key element for the clas- 
sification of T Tauri stars (Krautter et al., 1997). If such a star 
is present, it should be similar to a weak emission line T Tauri 
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Si3.0 sun ^13.2 Si 3.-1 S13.5 513.6 



Fig. 6. The spectrum of HR 6000 (full thick line) is compared with that of HD 175640 (dashed line) in a region crowded with lines 
not predicted by the synthetic spectrum (thin line). The labels indicate the only identifications in this spectral range. 



star (WTT). WTT stars are defined as low-luminosity TTauri 
stars with strong X-ray emissions and with an optical counterpart 
showing pre-main-sequence characteristics, in particular strong 
Li I absorption at 6707 A and very weak or also absent emissions 
in Ha. Furthermore, the radial velocity of the WWT stars must 
be consistent with that of the members of the molecular cloud to 
which they belong (Bertout, 1989). 

Upper panel of Fig. 8 shows that the observed spectrum 
is much better reproduced by a convolution of the synthetic 
spectrum of HR 6000 with a synthetic spectrum computed for 
Teff =3500 K, log g=4.Q and solar abundances. The computed 
spectrum of the cool star was not shifted in radial velocity rel- 
atively to that of HR 6000. In this case the contribution of Li i 
6707.7 improves the agreement between the computed spectrum 
and the observed spectrum whose continuum was lowered by 
1.25%. We note that the main contribution to the whole absorp- 
tion comes from the TiO lines whose labels are not plotted in 
Fig. 8, due to their huge number. All the considerations relative 
to Fig. 8 can be extended to the whole spectrum, starting from 
about 5800 A. 

To model the combined spectrum shown in the upper panel 
of Fig. 8 we assumed a luminosity ratio Li,/LTTauri='2-2Q, which 
is among the values suggested by van den Ancker et al. (1996). 
When the distance d=140+20pc adopted for HR6000 by van 
den Ancker et al. (1996) is replaced by the Hipparcos distance 
d=240±48pc, the stellar luminosity computed with formula (1) 
from van den Ancker et al. (1996) is L*=324Lo. In this case 
the luminosity of the T Tauri star becomes L7-7'a„ri=1.5Lo. For 
the computation of the binary spectrum the ratio of the stellar 
radii i^TTaJ^HRmm))^ was set to be equal to 0.8 in accordance 



with the Stefan-Boltzmann law. A change of the luminosity ra- 
tio affects only the level of the computed continuum. Several 
experiments made with other choices of the parameters with Teff 
ranging from 3500 K to 4500 K and log g ranging from 4.0 to 1 .5 
have shown that the above results are not significantly modified. 

In any case, although we can not exclude in an absolute way 
that the spectrum of HR 6000 is contaminated by that of a WTT 
star, we can not exclude as well that we simply observe spectral 
noise rather than weak features of a secondary spectrum. In fact, 
the inspection of UVES spectra of other CP stars observed dur- 
ing the same runs has shown the presence of weak features very 
similar to those we just discussed here for HR 6000. 



8. Conclusions 

The present study of HR6000 has led to the following results: 
the measured radial velocity indicates that HR 6000 most likely 
belongs to the Lupus cloud so that its age is of the order of 10^ 
years as it was found for the Lupus 3 cloud (James et al., 2006). 
Interstellar lines from the cloud can be observed in the spectrum 
from the ultraviolet to the visible, with displacements of the or- 
der of -0.05 -0.15 A from the stellar Unes. 

UVES spectra observed with a six months time interval do 
not show spectrum variability for both stellar and interstellar 
lines, but a low variability of 1.2 km s ' in the radial velocity 
and marginal variabilities of the weak spectral features which 
are however at the level of the noise. 

Different methods for the parameter determination have led 
to large differences in T^ff. While the gravity values are found 
in the range between 4.1 and 4.4 dex with errors of the order 
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178" 
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Modined 


predicted(K03)'' 


5130.0 


Fen 


5130.011 


5130.013 


-5.590 


-5.358 


-1.790 


5130.006 


1.172 


5130.18 


Fen 
Fell 


5130.177 
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5130.196 

CI OA 1 no 

5130.198 


-0.673 
-0.573 


5130.3 
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-6.352 




-1.252 






5130.45 


Fell 












5130.427 


-0.568 


5130.60 


Fen 
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5130.618 
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-0.312 
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-0.570 
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5131.706 
5131.788 


-5.127 
-5.023 




-0.127 
-0.723 


5131.878 


-0.727 


5132.0 


Fen 




5132.002 


-4.648 


-4.590 


-1.648 


5131.953 


-0.222 


5132.55 


Fen 












5132.554 


-0.190 


5133.59 


Fen 




5133.589 


-3.961 




-1.061 


5133.475 


0.220 


5134.07 


Fen 




5134.082 


-5.805 




-0.105 






5134.2 


? 
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5135.25 


Fen 
Fen 


5135.268 
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5135.187 


-6.671 




-2.671 








Fen 




5135.297 


-5.119 


-4.394 


-1.719 







5135.35 ? 

5135.95 ? 



''J78=Johansson (1978) 
*K03=Kurucz (2003) 
'^RU=Raassen & Uylings (1998) 



of 0.1-;-0.2dex, T^n changes from 1 3900+300 K, (if obtained 
from the UBVjS and uvby/? photometry), down to 12850±50K 
when determined from the Bahner profiles, and can become as 
low as 11200±300K when derived from the VIeJ photometry. 
Furthermore, as far as the Balmer profiles are concerned, the pa- 
rameters which well reproduce Hp, Hy, and predict too nar- 
row H„ wings. All these discrepancies in the parameter deter- 
minations could be accounted for by the presence of the Lupus 
cloud with could affect the colors in such a way to invalidate 
the standard reddening relations. Also a spectral contamination 
by a weak-emission T Tauri (WTT) companion, as suggested by 
van den Ancker et al. (1996), can not be excluded in an abso- 
lute way, although the only spectroscopic sign of its possible 
presence is a weak variable broad absorption at the position of 
Li 1 6707 A and a somewhat better agreement between observed 
and computed spectra longward of 6000 A for a combined syn- 
thetic spectrum obtained from that HR 6000 and one computed 
for retf=3500K, log ^=4.0, fM/H]=0.0. This companion could 
be either physically related with HR6000 to form a close bi- 
nary system or could be located in the foreground of HR 6000. 
Finally, the different parameters derived from the different de- 
terminations could be due to the extremely peculiar nature of 
HR 6000 which may make the classical model atmospheres in- 
adequate to reproduce all the observed stellar quantities. In fact, 
the abundance stratification inferred for some elements, heUum 
in particular, would require more refined models in which the 
hypothesis of constant abundances throughout the atmosphere is 
dropped. 

By using an ATLAS 12 model with parameters reff= 12850 K, 
log g=4.1,^=0.0kms"' we have computed a synthetic spectrum 
for HR 6000 from 3050 A to 9460 A and have compared it with 
the observed spectrum. This comparison has shown that most 
of the elements lighter than Ca are significantly underabundant, 
except for Be, Na, and P. The Si underabundance ([-2.9]) is re- 
markable because it is even lower than that of 46 Aql ([- 1 .0]) 



which was claimed by Sadakane et al. (2001) to be the lowest 
one found in HgMn stars. 

A striking peculiarity of HR 6000 is the lack of any over- 
abundance for heavy elements with Z>40, except for Xe and 
Hg. The line spectrum of Xeii is similar to that we observed 
in a prehminary analysis of UVES spectra of 46 Aql, another 
Xe Rich star: while most of the lines lie at the laboratory wave- 
length, some other Xe ii lines seem to be shifted by about -0.1 A 
from the predicted position. The same behaviour can be ob- 
served also in HD 175640, although to a less extent, owing to 
its lower Xe overabundance. Unfortunately, the small number of 
known transition probabilities and the ignorance of the Xe ii iso- 
topic structure due to the lack of atomic data put strong limita- 
tions in the study of this element in the spectra of the CP stars. 

The large overabundance of [-1-0.7] for iron generates a very 
rich fine spectrum in which numerous, still unclassified Fen 
Unes are observed. A very large number of lines probably due 
to Fe n remain unidentified. The similar iron overabundance of 
46 Aql ([+0.65]) gives rise to an impressive close resemblance 
between the spectra of the two stars. 

Abundances in the ultraviolet obtained from a re-analysis 
of the lUE spectra studied by Castelli et al. (1985) agree with 
abundances derived from the UVES spectra for all the elements, 
except for carbon and phosphorus. While the carbon identifica- 
tion in the optical spectrum is rather questionable owing to the 
weakness of the observed lines, the difference in the phosphorus 
abundance is similar to that yielded by the individual P n and P m 
lines. 

There are several signs of vertical abundance stratification in 
HR 6000 for He, P, Mn, and Fe. The peculiar shape of the He i 
profiles is discussed for the first time in this paper. The profiles 
have cores too intense as compared to the wings, a fact indicating 
strong vertical He abundance stratification (Dworetsky, 2004; 
Bohlender, 2005). Furthermore, the inferred He abundance de- 
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Fig. 7. The observed spectrum of HR 6000 (thick line) is compared with synthetic spectra (thin line) computed with the different 
choices of Fe ii log g/'s. The lower panel shows the synthetic spectrum as predicted by the Kurucz (2003) line lists when only 
transitions from observed levels are considered. Middle panel shows the synthetic spectrum when modified log gf's for several Fe ii 
lines are adopted (Table 6, col. 7). Upper panel displays the synthetic spectrum computed by considering Fe n transitions from both 
observed and predicted levels. 



creases with increasing wavelength ranging from an underabun- 
dance of [-0.8] at 4000 A to [-1.6] at 6000 A. 

For Mnii, the abundance shortward of the Balmer disconti- 
nuity is larger by about 0.6 dex than that longward of the Balmer 
discontinuity; for phosphorus, the Pii abundance is larger by 
0.2 dex than the P m abundance, but this discrepancy lies within 
the mean square error of the average abundances; for iron, the 
abundance from Fe ii high excitation lines is generally higher by 
about 0.2 dex than that from Fe ii low excitation lines. 

Similar to other studied HgMn stars, also HR6000 shows 
emission lines. The most numerous ones belong to Mnii, fol- 
lowed by those of Fe ii and Cr n. While the first two elements 
are overabundant in HR 6000, chromium has solar abundance. 
The star exhibits also isotopic anomalies for Hg and Ca. In both 
cases the most heavy isotope is the predominant one. 
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